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ABSTRACT 



We present photometric, astrometric and spectroscopic observations of the 
nearby (9.8 pc) low-mass binary G1569Bab (in turn being a companion to the 
early-M star G1569A), made with the Keck adaptive optics facility. Having 
observed G1569Bab since August 1999, we are able to see orbital motion and to 
determine the orbital parameters of the pair. We find the orbital period to be 
892 ± 25 days, the semi-major axis to be 0.90 ± 0.02 AU, the eccentricity to be 
0.32 ± 0.02 and the inclination of the system to be 34 ± 3 degrees (1-(t). The 
total mass is found to be O.123]'^o;q22^0 In addition, we have obtained 

low resolution (i? = 1500-1700) near-infrared spectra of each of the components 
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in the J- and K-bands. We determine the spectral types of the objects to be 
M8.5V (G1569Ba) and M9V (G1569Bb) with an uncertainty of half a subclass. 
We also present new J- and K-band photometry which allows us to accurately 
place the objects in the HR diagram. Most likely the binary system is comprised 
of two brown dwarfs with a mass ratio of 0.89 and with an age of approximately 
300 Myr. 

Subject headings: stars: low-mass, brown dwarfs — binaries: close — Hertzsprung- 
Russell diagram — astrometry — instrumentation: adaptive optics 

1. Introduction 

Brown dwarfs (BDs), despite sometimes being labeled "failed stars", are very interesting 
objects. According to our current understanding, BDs may represent the extreme low- 
mass end of star formation in which the mass is too small to sustain thermonuclear fusion. 
This stellar-substellar transition (also defined as the minimum mass at which the internal 
energy provided by nuclear burning quickly balances the gravitational contraction energy) is 
expected to occur around 0.075 (Baraffe et al. 1998) for objects with solar metallicities. 
This value is slightly lower in recent models, which predict a transition mass of 0.072 A/q 
(Chabrier & Baraffe 2000). Objects with masses below this limit never reach the stable 
hydrogen burning main sequence, but instead cool down as they age (Burrows et al. 1997; 
Chabrier et al. 2000) so that their surface temperatures and luminosities strongly depend on 
age as well as mass. As these cooling curves rely on poorly tested theoretical models, it is 
highly desirable to calibrate them with direct measurements. 

Obtaining dynamical masses for very low-mass (VLM, Af < 0.2 Ada-) stars and BDs 
from binaries is a challenging prospect, as so far there are no known eclipsing binary VLM 
stars or BDs. However, there are a number of known wide, non-eclipsing VLM and BD 
binaries (Kirkpatrick et al. 2001; Gizis, Kirkpatrick & Wilson 2001; Leinert et al. 2001; 
see also Reid et al. 2001a and references therein) that have been observed with a range of 
instruments and techniques, including ground-based infrared and optical imaging, speckle 
interferometry, adaptive optics (AO), and the Hubble Space Telescope. These binary systems 
promise to yield highly accurate dynamical masses, although they tend to have long periods 
and hence will require patience. 

In the absence of direct mass measurements for VLM/BD objects, one has to rely on 
indirect methods that may constrain mass ranges but that do not provide high-precision 
mass values. One very useful such technique is the lithium test (Magazzii, Martin & Rebolo 
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1993). Lithium is an element that is easily destroyed under the conditions prevalent in stellar 
interiors at temperatures slightly below those required for hydrogen burning. Objects more 
massive than ~ 0.060 Mq have their primordial lithium abundances depleted as long as they 
are fully convective. The lithium test has been used to confirm BD candidates in the Pleiades 
(e.g., Basri, Marcy & Graham 1996) and in the field (e.g., Martin, Basri & Zapatero Osorio 
1999; Kirkpatrick et al. 1999; Tinney Delfosse & Forveille 1997). 

G1569A is a nearby ((i = 9.8pc), chromospherically active late-type (M2.5V) star. For- 
rest, Skrutskie & Shure (1988) first reported a possible BD companion with a separation 
of 5" from the primary. Based on a low resolution spectrum Henry & Kirkpatrick (1990) 
classified this object as an M8.5 dwarf with a mass of 0.09 ± O.O2M0. More recently, Martm 
et al. (2000a) resolved the companion into two separate objects (G1569Ba and Bb) using 
AO observations with the Keck II telescope. They estimated the orbital period of the Ba- 
Bb binary to be around 3 years, and the total mass of the binary pair to be in the range 
0.09 — 0.15 Mq. Herein we present the results of extensive follow-up observations of this inter- 
esting pair, including improved photometry as well as near-lR spectroscopy and astrometry 
of the resolved binary components. We use the photometry to accurately place the objects 
in the Hertzsprung-Russell (HR) diagram, the spectroscopy to derive spectral types, and the 
astrometry to derive the total mass of the system from its orbital parameters. 



2. NIR photometry and astrometry 

2.1. The composite pair: G1569Bab 

Broad-band near- infrared photometry of the composite pair G1569Bab is available in 
the hterature (Forrest et al. 1988; Becklin & Zuckerman 1988). However, the measurement 
uncertainties claimed by the authors are too large for an accurate placement of these objects 
in the HR diagram or for direct comparison with theoretical evolutionary models. With the 
objective of improving the photometric data, we have collected J and Kghort direct images 
of the system G1569A and G1569Bab with the near-infrared camera (HgCdTe detector, 
256 X 256 elements) mounted at the Cassegrain focus of the 1.5-m Carlos Sanchez Telescope 
(GST, Teide Observatory) on 2001 February 8. The observations were performed through the 
"narrow-optics" of the instrument, which provides a pixel projection of 0.4" onto the sky. 
The atmospheric seeing conditions during the night of the observations were fairly stable 
around 1", which allowed us to easily separate the M2.5-type star from the pair G1569Bab. 
This latter object was not resolved into its two components. The total integration times 
were 5 s and 40 s in J and K short filters, respectively. A five-position dither pattern was 
used to obtain the images; each image consisted of 4 (J) or 8 {K) co-added exposures of 
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0.25 s (J) and 0.5 s (K) respectively. The dither pattern was repeated twice for the K-band 
observations. 

Dithered images were combined in order to obtain the sky background, which was later 
substracted from each single frame. Gl 569Bab is clearly detected in individual images, and 
we have obtained aperture photometry on each of them using PHOT in IRAF.^ Instru- 
mental magnitudes were placed on the UKIRT photometric system using observations of 
the standard star HD 136754 (Casali & Hawarden 1992), which was imaged with the same 
instrumental configuration just before and immediately after our target. Both the science 
target and the standard star were observed at similar air masses. The photometric error of 
the calibration was ±0.03 mag in both filters. Kghort displays a different bandpass compared 
to Kukirt] the transformation between these two filters is not well defined yet, albeit for 
objects as red as Gl 569Bab it has been estimated at Kshort-KuKiiar — 0.035 (Hodgkin et al. 
1999, and references therein). We have applied this correction to our photometry as well as 
the relations given in Leggett, Allard & Hauschildt (1998) to convert UKIRT data into the 
CIT photometric system. The final average magnitudes derived for G1569Bab are given in 
Table 1, where the photometric errors listed correspond to typical 1 cr uncertainties of single 
measurements. 

The astrometry of Gl 569Bab relative to the bright primary Gl 569A as measured on the 
GST data (MJD = 51948.202) is the following: angular separation of 4."890 ± 0."040, and 
position angle of 30° ±3°. We note that these values differ from those published in Forrest 
et al. (1988) by more than 2(7, providing evidence of the orbital motion of G1569Bab around 
the M2.5-type star in the time interval of roughly 15 years. 

2.2. Adaptive Optics Imaging of G1569Ba-Bb pair 

G1569Bab was observed on 7 occasions between 1999 August and 2001 May with the 
Keck 11 AO system (Wizinowich et al. 1988). The first 3 observations made use of the 
KGAM camera with a NIGMOS-3 infrared array. The later observations made use of the 
slit- viewing camera (SGAM) associated with the NIRSPEC instrument (McLean et al. 1998). 
SGAM uses a PIGNIG HgGdTe array. For both cameras the pixel scale was 0.0175" and the 
field of view was 4.48" x 4.48", except for the 2001 May observation when the SGAM pixel 
scale was changed to 0.0168". Exposures were generally obtained in the X'-band, except 
for 2000 Febuary when the observations were taken in the J-band. Exact exposure times 
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varied, but typically consisted of 30 coadded 2-second exposures. We used a 1% transmission 
neutral density filter in the beam to prevent saturating the bright primary star. Flat-field 
correction was performed using twilight exposures, while sky subtraction made use of images 
of an adjacent field observed immediately after G1569B. Corrected seeing varied between 
0.05" and 0.08". For three of the observations the primary star (G1569A) was also in the 
field of view, providing an in-field astrometric reference. For the other observations the 
primary was either not observed or was saturated. No photometric standard stars were 
observed in any of the epochs, so we cannot provide absolute photometric calibrations for 
the AO observations. Figure 1 shows the resulting images of the pair G1569Bab at six 
different epochs. This figure clearly demonstrates that this system is resolved into a binary 
and that orbital motion is evident. 

We used the DAOPHOT package (in IRAF) for data reduction and analysis. The 
point spread function of the objects in each frame were fitted with an elliptical Gaussian 
function, providing relative astrometry (Table 2) and photometry. The relative photometry 
of G1569Ba and G1569Bb was derived by computing the ratio of the amplitudes of the best 
fitting Gaussians (Table 1). Gl 569Bb is fainter by 0.51 ± 0.02 mag and 0.41 ± 0.03 mag in the 
J- and X-bands, respectively. This makes this object redder in (J — K) by 0.10 ± 0.04 mag. 
The error bars take into account the dispersion observed from image to image, and from one 
observing run to another. We do not find a significant relative photometric variability in 
any member of the pair within 3 cr the uncertainties. With the relative brightness of the two 
components and the combined fiux known, it is possible to derive the individual absolute 
magnitudes of Gl 569Ba and Gl 569Bb. We fist in Table 1 the resulting decomposition for the 
J- and X-bands. The corresponding error bars incorporate the photometric uncertainties of 
the combined system Gl 569Bab and the uncertainties of the relative photometry. We are 
confident that the latter is determined with a higher accuracy. 



3. Low resolution NIR spectra 

We have obtained low resolution spectra of G1569Ba and G1569Bb in the J- (1.158- 
1.368 /im) and i^-bands (1.992-2.420 /im) using the cross-dispersion spectrograph NIRSPEC 
and the AO facility at the Keck 11 telescope. The data were collected on 2000 June 20. The 
raw seeing and transparency conditions were very good during the observations, and the AO 
correction applied to the primary star G1569A provided well resolved images of the binary 
system G1569Bab (AO corrected seeing of 0.05"). NIRSPEC in its spectroscopic mode is 
equipped with an Aladdin InSb 1024 x 1024 detector with a pixel projecting 0.0185" onto 
the sky. For the present study, we selected the low resolution spectroscopic mode which 
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provides nominal dispersions of 2.8 A/pix and 4.2 A/pix in the J- and X-bands, respectively. 
The 3 pixel-wide slit was aligned with the two components of the binary (PA~139deg) so 
that both targets were observed simultaneously. 

Total exposure times were 240 sec and 400 sec for the J-band and X-band spectra, 
respectively. The observing strategy employed was as follows: 6 (J) and 10 (K) individual 
integrations of 20 sec (J) and 10 sec (K) each at two different positions along the entrance 
slit separated by about 1.8". This procedure was repeated twice in the K-hand. In order 
to remove telluric absorptions due to the Earth's atmosphere, the near-infrared featureless 
AOV-typc star HR5567 was observed very close in time and in air mass (within 0.05 air 
masses). Calibration images (argon arc lamp emission spectra and white- light spectra) were 
systematically taken after observing each source. 

Raw data were reduced following conventional techniques in the near- infrared. Nodded 
images were subtracted to remove the sky background and dark current. The spectra of the 
sources and of the cahbration lamps were then extracted using subroutines of the TWOD- 
SPEC package available in IRAF. The extraction apertures of G1569Ba and G1569Bb were 
selected so that cross-contamination was less than 10%. The extracted spectra of the sources 
were divided by their corresponding normalized extracted flat-fields, and calibrated in wave- 
length. The la dispersion of the fourth-order polynomial fit was 0.4 A and 1.0 A in the J 
and K spectra, respectively. The hydrogen Pf3 absorption line at 1.2818 /xm and the B7 
absorption line at 2.1655 //m in the spectra of HR5567 were interpolated before they were 
used for division into the corresponding science spectra. We are confident that the science 
spectra have good cancellation of atmospheric features. To complete the data reduction, 
we multiplied the spectra of our targets by the black body spectrum for the temperature of 
9480 K, which corresponds to the AOV class (Allen 2000). 



3.1. Spectral types, atomic and molecular features 

The resultant average spectra with a resolution of i?=1500 in J and i?=1700 in K 
are depicted in Fig. 2. The strongest molecular and atomic features are indicated following 
the identifications provided by Jones et al. (1996) and McLean et al. (2000). The spectra of 
both components, Gl 569Ba and Gl 569Bb, are indeed very similar. The composite spectrum 
of G1569Bab in the optical has been previously studied by Henry & Kirkpatrick (1990) 
and Kirkpatrick, Henry & McCarthy (1991), who derived a dwarf spectral type of M8.5. 
In addition, this object is listed in the Table 1 of Kirkpatrick et al. (1991) as a primary 
dwarf spectral standard. Our data agree with this measurement for the bright component 
G1569Ba, and also provide evidence that G1569Bb is not significantly cooler. This is fully 
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consistent with the photometry presented above. 

We have obtained equivalent widths of the strongest observed atomic absorptions of K I 
and Na I in the spectra; the measurements are given in Table 3. Due to the low resolution of 
our data, the majority of these lines are considerably blended with other spectral features, 
e.g. the Kl line at 1.2435 //m is contaminated by a strong molecular band of FeH. The 
values in Tabic 3 have been extracted adopting the base of each line as the continuum. 
We find typical standard deviations in equivalent width close to 10% over the reasonable 
range of possible continua. Although this procedure does not give an absolute equivalent 
width, it is commonly used by different authors, and allows us to compare our values with 
those published in the literature. We have also measured the strengths of the H2O band 
at 1.330 //m and the CO band at 2.294 //m in a way similar to that described in McLean 
et al. (2000), Reid et al. (2001b) and Jones et al. (1994). Our measurements are listed in 
Table 3 with uncertainties of about 5%. All these values are comparable to those obtained 
from similar spectral type field stars, which suggests that neither G1569Ba nor G1569Bb 
have very discrepant metallicities or gravity. 

We note that the equivalent widths of the K I hues and the H2O and CO absorptions 
in G1569Bb appear to be shghtly larger than in G1569Ba, while the Nal hues in the K- 
band spectrum are smaller. This trend is observed for decreasing temperatures (Jones et al. 
1994), and clearly indicates the cooler nature of G1569Bb. By fitting a polynomial spectral 
type-equivalent width relation to the data available for spectral standard stars (see Reid 
et al. 2001b), we conclude that the differences between Ba & Bb in our measurements are 
consistent with G1569Bb being half a subclass cooler. This would make G1569Bb an M9- 
dwarf (±0.5 subclasses). 



3.2. Radial velocities 

We used our low resolution near-lR spectra taken on 2000 June 20 (MJD = 51715.365) to 
compute the relative radial velocity of Gl 569Bb and Gl 569Ba via Fourier cross-correlation. 
Because the spectroscopic data have been corrected for telluric lines, we do not expect these 
lines to be a large source of uncertainty. Unfortunately, no spectra were taken of the primary 
star G1569A, so we cannot determine the relative radial velocity of the pair with respect it. 
The velocity dispersion of the data is rather poor (lpixel~66kms~^ in J and ~57kms~^ 
in K). Nevertheless, the cross-correlation technique was able to achieve precisions of about 
1/4 pixel, so we obtained a relative radial velocity accurate to about 15kms^^. We verified 
this by cross-correlating individual spectra of each component against itself. The relative 
velocity (G1569Bb cross-correlated with G1569Ba) we measure is 25kms~^ in J and 6kms~^ 
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in with an average value of 15.5 km s~^. The peak-to-peak radial velocity variation of 
the system on the basis of the orbital solution presented in next section is around 14kms~^; 
our measurement is consistent within the error bar with the expected value at the epoch 
of the observations. However, this error bar is rather large and prevents us from making 
further analysis (like the presence of invisible companions). The maximum peak of the cross- 
correlation function is in the range 0.93-0.97, which indicates the similarity and the high 
signal-to-noise ratio of the spectra. 

4. Orbit determination and total mass of the pair 

We determined the apparent orbit of the binary pair G1569Bab by fitting a Keplerian 
model to the relative astrometric data shown in Table 2. As the fit is non-linear in the 
orbital elements, we made use of a gradient-following fitting routine (Press 1992) to find the 
optimal (in a chi-squared, x^, sense) orbital parameter set. As there may be many local 
minima in the manifold, and our gradient following routine moves strictly downhill, we 
ensured that it found the global minimum by starting it at a range of different locations 
in parameter space. The best-fit parameters are given in Table 4, and the visual orbit is 
shown together with the astrometry in Fig. 3. The rms of the residuals is 0.0024", and the 
residuals do not show any long-term drifts. Uncertainties (assuming normal errors) were 
estimated from the covariance matrix of the fit, scaled by the reduced chi-squared (0.96). 
However, as apparent in Fig. 4 the best-fit values for the semi-major axis and period are 
correlated, leading to a lower uncertainty in the total mass. The mass uncertainty was 
estimated by finding contours in parameter space where the x^ was increased by 2.3, 6.2 
and 11.8 respectively, corresponding to the lo", 2(T and 3(7 contours. The Icr, 2(7 and 3(7 
mass ranges are 0.114-0.135, 0.107-0.142 and 0.101-0.150 Mq respectively, while the best-fit 
values for the period and semi-major axis correspond to a total mass of 0.123 Mq. From the 
relative photometry we infer that the mass ratio of the binary is close to, but not exactly, 
equal. Hence the 2a upper mass limit of the secondary G1569Bb is less than 0.071 Mq, i.e. 
below the hydrogen burning mass threshold, making it a likely brown dwarf. 

It is important to note that although we were able to obtain reliable relative astrometry 
of the Ba-Bb pair (over a separation of ~ 0.1"), the uncertainties in plate scale and orientation 
and the saturation of the bright star in some images were such that we were unable to reliably 
measure the orbital motion of each component of the Ba-Bb pair with respect to a separate 
reference, i.e. G1569A (located ~ 5" away). Hence, while we can determine the relative 
orbit of the Ba-Bb pair to a high degree of precision, we cannot astrometrically determine 
the mass ratio of the two components. Further studies are needed to confirm or discard the 
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substellar nature of the primary Gl 569Ba. We will combine our photometry (absolute and 
relative), the total mass of the system and additional information available in the literature 
to compare the binary with the most recent evolutionary models. 

5. Discussion 

5.1. Color- magnitude diagram 

Figure 5 depicts the location of the Gl 569Bab pair in the near-infrared color-magnitude 
diagram. To convert their observed magnitudes into absolute magnitudes wc have used the 
astrometric parallax provided by Hipparcos (0. "10191 ± 0. "00167, Ferryman et al. 1997), 
which is very similar to previous astrometric measurements (Heintz 1991). Also shown in 
this figure are the locations of very late-type dwarfs in the Pleiades cluster (~120 Myr, Basri 
et al. 1996; Martin et al. 1998; Stauffer, Schultz & Kirkpatrick 1998, we use a distance of 
120 pc) which have photometry available in the hterature (Festin 1998; Martin et al. 2000b), 
and of objects in the field. Absolute magnitudes and colors of M-typc field standard stars 
have been taken from tables published in Kirkpatrick & McGarthy (1994) and Leggett et al. 
(1998). For L-type field dwarfs wc have adopted the average near-infrared colors provided in 
Kirkpatrick et al. (2000), and have averaged absolute K magnitudes for those objects with 
parallax available in the literature (see Reid et al. 2000, 2001b; Kirkpatrick et al. 2000). 

Overplotted onto the observed data in Fig. 5 are the 0.5, 1.0 and 5.0 Gyr theoretical 
solar composition isochrones from the evolutionary models of the Lyon group (Chabrier et 
al. 2000; BarafTc et al. 1998) and from those of the Arizona group (Burrows et al. 1997). 
We have adopted solar metallicity in our studies because the photospheric abundance of the 
bright star G1569A has been determined to be very close to solar ([Fe/H] =-0.15, Zboril & 
Byrne 1998). Although the Lyon models do provide magnitudes and colors in the filters of 
interest, we preferred to compute them from the predicted luminosity and effective temper- 
ature at a given mass and age. This allows a direct comparison of the two sets of interior 
models and minimizes the effects of possible errors in the model atmosphere synthesis. We 
converted effective temperatures (Tgfj) into colors by using the temperature scales of Leggett 
et al. (1998) for M-dwarfs and of Basri et al. (2000) for late-M and L-dwarfs. These two 
temperature scales should be consistent and complementary with each other as the authors 
make use of the same atmosphere models to obtain their results. We derived absolute K 
magnitudes from theoretical luminosities by using the bolometric correction as a function 
of spectral type (i.e. color, T^s) given in Leggett et al. (2000) down to mid-M classes, and 
in Reid et al. (2001b) for cooler types. To summarize, the second order polynomial fits 
(1700 < Teff < 3500 K) we used are as follows: 
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(J - X) = 6.423 - 3.49 x lO-^^Teg + 5.41 x lO-^s rms = 0.04 mag 

BCk = 5.745 - 1.46 x 10-=^Teff + 1.57 x lO-^ff = 0.06 mag 

Isochrones in Fig. 5 are plotted for T^s < 2900 K, which roughly corresponds to masses smaller 
than 0.1 at ages around 1 Gyr. 

From Fig. 5 we can see that the evolutionary models nicely reproduce the trend delin- 
eated by field objects, except for the reddest colors {J — K > 1.5) where models apparently 
predict brighter magnitudes. Of the two sets of isochrones, the Lyon 1-5 Gyr models seem 
to produce a better fit to the observed data in the field. The difference in color between 
Gl 569Ba and Gl 569Bb is consistent with the spectral types of the objects. Within 1 a the 
uncertainties of our JK photometry, the location of the pair is well matched by isochrones 
in the age interval 0.2-1.0 Gyr. This indicates a young age for the multiple system, a result 
which is compatible with the elevated X-ray emission of the "single" M2.5-type primary 
(Pallavicini, Tagliaferri & Stella 1990; Huensch et al. 1999), with the system belonging to 
the young galactic disk as inferred from its kinematics (Reid, Hawley & Gizis 1995), with 
the large rotation rate measured for the star G1569A (Marcy & Chen 1992), as well as with 
the late-M spectral type-lithium-age relationships (Magazzii et al. 1993; see Bildsten et al. 
1997). 

Lithium is detected in M8-M9 Pleiades BDs (Rebolo et al. 1996; Stauffcr et al. 1998), 
whereas older and slightly more massive objects have depleted it very efficiently. Thus, 
lithium non-detections in very late-M type objects necessarily imply ages older than the 
Pleiades. No hthium feature is observed in the composite optical spectrum of G1569Bab 
(Magazzii et al. 1993), thus implying that the binary is older than 0.12 Gyr. The it!-band 
spectroscopic data shown in Magazzii et al. (1993) have poor signal-to-noise ratio, and are 
dominated by the bright and more massive component as it contributes twice as much flux 
as does the fainter component. We will discuss later how the age of the system can be 
constrained to a much smaller range using results presented in this paper. 

We do not flnd from our near-infrared photometry strong evidence for the possible 

binary nature of G1569Ba as claimed by Martin et al. (2000a) on the basis of their H vs 
H — K color-magnitude diagram as well as by Kcnworthy et al . (2001) on the basis of their 
J — K colors. The relative position of the two components of the pair in Fig. 5 reasonably fits 
the location of field dwarfs even within 1 a the error bars. If G1569Ba is a binary itself, the 
smaller companion has to be at least a factor five less luminous. We have combined our AO 
K-hand images to look for any possible companion. We place a 3(j limit at K — 16.5ma.g 
(0.015-0.02 Mq) on the brightness of a possible companion at distances greater than 0.25", 
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and less than 2"- half the size of the AO detector. We cannot discard, however, the presence 
of extremely faint and less massive objects around any of the components which our AO 
observations have not been able to detect /resolve. Follow-up high resolution spectroscopy 
and/or very detailed analysis of the orbital motion of the pair may reveal the presence of 
close-in giant planets. 

5.2. The HR diagram and substelleirity 

Figures 6 and 7 show the location of Gl 569Ba and Gl 569Bb in the HR diagram (lumi- 
nosity as a function of effective temperature) and provide a comparison with state-of-the-art 
evolutionary models by the Lyon group (Ghabrier et al. 2000) and by the Arizona group 
(Burrows et al. 1997). Solar- mctallicity abundance isochrones of ages 120, 300, 500 Myr and 
1 Gyr, and evolutionary tracks of masses in the interval 0.04-0.09 are shown in these 
figures. We use the most recent determination of the substellar mass limit at 0.072 Mq to 
define the stellar-substellar borderline. Because this value is the smallest of those available 
in the hterature (see e.g. Grossman, Hays & Graboske 1974; D'Antona & Mazzitelli 1994; 
Burrows et al. 1993), our conclusions on substellarity will be conservative. We indicate in the 
figures the substellar mass boundary and the location of the 50% depletion limit of lithium 
burning predicted by the two sets of models. We have obtained the luminosity and effective 
temperatures of our targets as explained above. The third-order polynomial fit that gives 
temperatures as a function of the observed (J — K) color is the following: 

Teff = 7744.6 - 9488.4( J - K) + 5509.9( J - Kf - 1135.7( J - Kf rms = 50 K 

This fit has been calculated for colors in the range 0.85-2.06 (spectral types M6-L6), and 
is based on the temperature cahbrations provided by Basri et al. (2000) and Leggett et 
al. (2000). Bolometric corrections in the J- and X-bands from Rcid et al. (2001b) and 
Bessell, Gastelli & Plcz (1998), respectively, have been used to transform magnitudes into 
bolometric luminosities. The values we derive for the pair arc listed in Table 1. The error 
bars in luminosity take into account the uncertainty of the distance modulus (Hipparcos) 
and the photometric uncertainties, leading to a total uncertainty of ±0.07 dex. The error 
bars assigned to the effective temperatures come from the uncertainty in the colors alone. 
These are in general a factor 2 larger than the rms of the polynomial fit describing the 
temperature calibration. 

Models should be able to provide explanations to all physical properties so far known 
for G1569Ba and G1569Bb, i.e., photometry, the total mass of the pair and the destruction 
of lithium. From Fig. 6 we observe that the location of Gl 569Ba is consistent with severe 
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lithium depletion, in agreement with available optical spectroscopic observations. According 
to the Lyon models, even the fainter companion G1569Bb has depleted its lithium. The 
likely age of the system is in the range 0.2-1 Gyr, but only with younger ages is it possible to 
reproduce the astrometric total mass derived for the pair. Therefore, the real constraint to 
the age of the system is given by the total mass rather than by the error bars in the figure. 
For the age of 300 Myr, the binary would be formed by objects of 0.069 (G1569Ba) and 
0.059 Mq (G1569Bb), in good agreement with the mean orbital solution. Both masses are 
below the stellar-substellar borderline and thus the two components would be brown dwarfs. 
For shghtly older ages, individual masses would become larger, as would the total mass. On 
the basis of the largest possible astrometric total mass value at the 3 a level, the pair would 
be made up of an object on the substellar borderline with 0.078 and a 0.070 M0-BD 
at the age of 500 Myr. We tabulate the possible masses of the pair as a function of age in 
Table 5; these estimations rely on the Lyon models. 

The comparison with the Arizona models shown in Fig. 7 also yields very young ages 
(120-500 Myr) for the system, and thus also small masses for each of the components. To 
this point, theory and some observations seem to be consistent. However, the Arizona 
evolutionary models predict hotter temperatures (by about 100 K) around the substellar 
mass limit at a given age than do the Lyon models, and hence according to these models 
both Gl 569Ba and Gl 569Bb should have preserved a considerable amount of lithium in their 
atmospheres. Only by assuming the highest temperatures allowed by the error bars can the 
Arizona models account for the observed lithium depletion in Gl 569Bab. This would move 
the pair to a location between 300 Myr and 500 Myr in Fig. 7. In order to be consistent with 
the additional restriction of the astrometric total mass, ages in the interval 300-500 Myr are 
required; the resulting masses are 0.055-0.075 for G1569Ba and of 0.048-0.068 M0 for 
G1569Bb. 

Collecting evidence from the orbital solution, photometry, spectroscopy, and the com- 
parison with evolutionary models, the most likely scenario of the binary Gl 569Bab is: two 
very late M-type BDs with masses of 0.055-0.078 M© and 0.048-0.070 M© in a close orbit, 
in turn orbiting an early M-type 0.5 Mq star, the whole system with an age in the range 
250-500 Myr. Such young ages found in nearby objects are not surprising since it now seems 
that the Sun is located close to a region that was the site of substantial amounts of recent 
stellar formation (Zuckerman & Webb 2000). 

Figure 8 portrays the mass-luminosity relationship for different ages as given by the 
evolutionary models. The two members of the pair are plotted with error bars indicating 
the uncertainty in luminosity and the hkely mass range of each component. Masses that 
have depleted hthium by a factor 2 are also incorporated into the figure. According to the 
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Arizona models, G1569Bb may have preserved lithium in its atmosphere, whereas this is 
quite unlikely based on the Lyon models. Therefore, lithium observations of this BD are 
needed in order to discriminate which model reproduces neatly the properties of the pair. In 
addition, precise radial velocity measurements of each component will lead to an accurate 
determination of the individual masses, and thus will also constitute a better constraint on 
the models. Nevertheless, models do not appear to be far from reproducing the observational 
properties of Gl 569Bab. Delfosse et al. (2000) show that the mass-luminosity relationship 
given by the Lyon models reasonably describes the low-mass stellar regime in the field. The 
pair Gl 569Bab has lower masses that belong to the substellar regime. Brown dwarfs around 
stars have been reported in the recent years (see Table 5 in Reid et al. 2001a for a compilation 
of the complete list), but to our knowledge none of them has been proved to be a binary 
itself. Gl 569Bab turns out to be the first confirmed resolved binary BD as a companion to 
a star. 

The distance to G1569 implies a physical separation of 49 AU between the M2.5-type 
star and the substellar pair. This large separation and the high mass ratio 0.135 and 
0.120) between the star and each of the BDs favors the fragmentation of a self-gravitational 
collapsing molecular cloud as the most plausible explanation for the formation of the system 
(Boss 2000; Bodenheimer 1998). Whether each component of the pair G1569Bab originated 
from a second fragmentation and collapse process of a small cloud core is not clear (the 
physical separation is 0.92 AU, and the mass ratio of the pair is ^"^0.89). The activity 
of the nascent low-mass star when it was gaining mass and becoming more luminous may 
have caused the disruption of the less massive collapsing core into two close substellar objects 
before the hydrostatic core could build up enough mass to eventually start hydrogen burning. 
Energetic outflows and jets up to thousands of AU in length have been detected in low-mass 
stars of very young star forming regions (e.g., Reipurth et al. 2000; Fridlund & Liseau 1998). 
We cannot discard the possibility, however, that the protoplanctary disk around the star 
might have also played an important role in the origin of the companions. Disks extending 
up to several hundred AU are known to exist around stars (Bruhweiler et al. 1997). Clearly, 
finding other similar systems will, in addition to providing additional dynamical masses, also 
contribute to our knowledge of the genesis of such interesting multiple low-mass binaries. 



6. Conclusions 

We have obtained new observations of the Gl 569Bab pair (low-mass binary companion 
at a wide separation from the M2.5-type star G1569A), which have allowed us to derive 
the spectral types of each component and the orbital parameters of the system. We find 
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that the total mass of the low-mass binary is O.1231;q;q22M0 (3-(t) with two detected compo- 
nents of M8.5 and M9 spectral types (half a subclass uncertainty) completing one eccentric 
(e = 0.32 ± 0.02) orbit every 892 ±25 days. We have also acquired new J and K near-lR 
photometry in order to locate Gl 569Ba and Gl 569Bb in the HR diagram and compare them 
with the most recent evolutionary models by the Lyon group (Chabrier et al. 2000) and the 
Arizona group (Burrows et al. 1997). The pair is likely formed by two solar metallicity young 
brown dwarfs with masses in the interval 0.055-0.078 Mq (G1569Ba) and 0.048-0.070 
(G1569Bb) at the young ages of 250-500 Myr. Our adaptive optics images taken with the 
Keck II telescope exclude the presence of any other resolved companion with K magnitudes 
brighter than 16.5 (3(t) at separations of 0.25" up to 2" from G1569Bab. This detection 
limit corresponds to masses around 0.015-0.02 Mq for the possible age range of the system. 
Further radial velocity and astrometric measurements will be very valuable to detect giant 
planets, as well as to provide individual masses for each of the members of the pair. 

We thank N. Reid for valuable comments and discussion. Data presented herein were 
obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among 
the California Institute of Technology, the University of California and the National Aeronau- 
tics and Space Administration (the Observatory was made possible by the generous financial 
support of the W. M. Keck Foundation); and at the Carlos Sanchez Telescope (GST) op- 
erated on the island of Tenerife in the Spanish Observatorio del Teide of the Instituto de 
Astrofisica de Canarias. This research has made use of the Simbad database, operated at 
CDS, Strasbourg, France. We are thankful to I. Baraffe and A. Burrows for providing us 
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Fig. 1. — Contour images of the Gl 569Ba-Bb pair (near-IR filter is given in brackets) showing 
the orbital motion. These data have been obtained with the Adaptive Optics facility of the 
Keck II telescope and with the KCAM (first 2 epochs) and SCAM/NIRSPEC (last 5 epochs) 
instruments. Data of 1999 August, when the binary was resolved for the first time, were 
presented in Martin et al. (2000a). 
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Fig. 2. — J-band (top panel) and i^-band (lower panel) NIRSPEC spectra of G1569Ba 
and Gl 569Bb obtained using the AO system of the Keck II telescope. Some features have 
been identified after Jones et al. (1994) and McLean et al. (2000). The spectra have been 
normalized to unity at 1.29 //m and at 2.19 jjtm. An offset has been added to G1569Ba's data 
for clarity. 
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Fig. 3. — The relative astrometry of the G1569Ba-Bb pair, together with the best-fit orbit 
(dotted ellipse). Error-bar crosses denote measurements and circles indicate the predicted 
location on the orbit at the time of the observations. North is up and East is to the left. 
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Fig. 4. — The ^ ^ function of period and semi-major axis. All other parameters are 
selected to provide the lowest x^- The contours give the la, 2a and 3a uncertainties in 
the two parameters. The two diagonal hues correspond to combinations of period and semi- 
major axis giving a total mass of 0.150 Mq and 0.101 Mq, respectively. These arc the upper 
and lower limits that we have adopted for the total mass of the binary Gl 569Bab. The cross 
indicates the preferred solution. 
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Fig. 5. — Infrared color-magnitude diagram displaying the location of G1569Bab (combined 
light) and each of the two components resolved with the Keck II Adaptive Optics system 
(filled circles). Pleiades members (~120Myr) are plotted with open circles, and the location 
delineated by field dwarfs with known parallax is shown with a thick full line. Isochrones 
provided by the Lyon group (0.5, 1.0 and 5.0 Gyr, Chabrier et al. 2000 — dashed lines) and 
by the Arizona group (0.5 Gyr, Burrows et al. 1997 — dotted line) are also overplotted in 
the diagram. Masses of 0.09 M© (crosses), 0.08 Af0 (plus-signs), 0.072 (open triangles), 
0.060 M0 (asterisks) and 0.055 M© (diamonds) arc marked with crosses on the isochrones. 
We indicate spectral types as a function of the {J — K) color on the bottom of the figure. 
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Fig. 6. — HR diagram illustrating the location of Gl 569Ba and Gl 569Bb (filled dots) in com- 
parison with theoretical evolutionary tracks of constant mass (dotted lines) and isochrones 
(full lines) from the Lyon group (Chabrier et al. 2000). The track corresponding to the 
substellar mass limit at 0.072 M0 is shown with a thicker dotted line. Masses in solar units 
are labelled on the upper part of the diagram, and ages for the isochrones are indicated to 
the right. The thick dashed line indicates the 50% lithium depletion limit predicted by the 
Lyon models. Objects to the left have severely depleted lithium, whereas objects to the right 
still preserve a significant amount of this element. Solar abundance has been assumed in 
generating this figure. 
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Fig. 7. — HR diagram illustrating the location of Gl 569Ba and Gl 569Bb (filled dots) in com- 
parison with theoretical evolutionary tracks of constant mass (dotted lines) and isochrones 
(full lines) from the Arizona group (Burrows et al. 1997). See the caption of Fig. 6 for further 
details. Here, the 50% lithium preservation line (thick dashed line) is taken from the Arizona 
models for consistency. 
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Fig. 8. — Mass-luminosity relationship for different ages (120, 300, 400 and 500 Myr, and 
1 Gyr) according to the Lyon models (Chabrier et al. 2000 — full hnes) and to the Arizona 
models (Burrows et al. 1997 — dashed lines). The stellar-substellar transition mass range at 
0.072-0.080 M0 is indicated by the vertical dotted hnes. The dot-dot-dot-dashed line (Lyon) 
and the dot-dashed line (Arizona) mark when 50% of the lithium has burned. Lithium 
preservation occurs to the right of these lines. The error bars assigned to G1569Ba and 
G1569Bb (filled circles) correspond to the uncertainty in luminosity, and to the likely mass 
range (S-cr) derived for each object. 
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Table 1: Photometry (CIT system) and physical parameters of G1569Ba and G1569Bb 



Object J K J-K log L/Lq ^eff (K) M/Mq 

G1569Bab 10.61 ±0.05 9.45 ±0.05 1.16±0.07 -3.17±0.07 0.101-0.150 

G1569Ba 11.14±0.07 10.02±0.08 1.12±0.10 -3.39±0.07 2440±100 0.055-0.078 

G1569Bb 11.65 ±0.07 10.43 ±0.08 1.22 ±0.10 -3.56 ±0.07 2305 ± 100 0.048-0.070 



Note. — The relative photometry between G1569Ba and G1569Bb is known to a better accuracy (see 
text). 



Table 2: Astrometry of G1569Ba-Bb 



Date 


Epoch 


Separation 


P.A. 




(MJD) 


(") 


n 


1999 Aug 29 


51419.270 


0.101 ±0.001 


46.8 ±3 


2000 Feb 18 


51592.600 


0.092 ±0.001 


98.2 ±3 


2000 Feb 25 


51599.644 


0.090 ±0.001 


100.4 ±2 


2000 Jun 20 


51715.408 


0.076 ±0.003 


138.6 ±2 


2000 Aug 24 


51780.283 


0.059 ± 0.001 


178.4 ± 2 


2001 Jan 09 


51918.665 


0.073 ± 0.002 


291.4 ±2 


2001 May 10 


52039.410 


0.097 ±0.001 


341.1 ±3 
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Table 3: Kl and Nal equivalent widths (A) and the strengths of the H2O and CO bands 







K I wavelength (//m) 


Na I wavelength (//m) 


HaO^ 


CO^ 


Object 


SpT 


1.169^ 1.177^ 1.244'= 1.253^= 


2.206^= 2.209^= 


1.330^= 


2.294^ 


G1569Ba 


M8.5 


6.5 7.5 4.7 4.8 


1.4 0.70 


0.72 


1.21 


G1569Bb 


M9.0 


6.5 7.7 5.2 5.1 


0.9 0.50 


0.70 


1.24 



"Ratio of the average flux in a 0.02 /im window centred on 1.34 /xm and on 1.29 y^m (Reid et al. 2001a). 
''Ratio of the average flux in a 0.06 /xm window centred on 2.25 iim and on 2.33 /xm (Jones et al. 1994). 
'^All wavelengths in /zm. 

Note. — Uncertainties are ±0.5 for the spectral classification, 10% for equivalent widths and 5% for the 
flux ratios. 



Table 4: Orbital parameters of Gl 569Ba-Bb 



Parameter 




Value 


Period, P 


892 


± 


25 day 


Eccentricity, e 


0.32 


± 


0.02 


Semi- Major Axis, a 


0.90 


± 


0.02 AU 


Inclination, i 


34 


± 


3 deg 


Arg. Periapsis, cu 


76 


± 


4 deg 


Long, of Ascending Node, D, 


141 


± 


4 deg 


Epoch (MJD), T 


51820 


± 


4 day 



Note. — Uncertainties are 1-a. Note that the toal mass uncertainty is smaller than this indicates, see 
Fig. 4 and text. 
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Table 5: Likely ages and individual masses of G1569Ba and G1569Bb 



Age 


G1569Ba 


G1569Bb 


Total mass 


(Myr) 


(Me) 


(Me) 


(Mo) 


200 


0.055 


0.048^^ 


0.103 


250 


0.061 


0.053^ 


0.114 


300'^ 


0.069^ 


0.059^ 


0.128"^ 


400 


0.072 


0.065 


0.137 


500 


0.078 


0.070 


0.148 


700^ 


0.082^ 


0.074^= 


O.ISG'^ 



"Should have (partially) preserved lithium. 

''These estimates match the mean astrometric orbital solution. 

'^The total mass inferred for this age is even beyond the 3 a uncertainty of our astrometric solution. 

Note. — Based on the Lyon models (Chabrier et al. 2000; Baraffe et al. 1998). The uncertainty of 
individual mass models is ±0.002 Mq. 



